Submitted to The Astro-physical Journal 



X-ray, Radio, and Optical Observations of the Putative Pulsar in the 

Supernova Remnant CTA 1 

J. P. Halpern, E. V. Gotthelf, F. Camilo, D. J. Helfand 

T^j- . Columbia Astrophysics Laboratory, Columbia University, 550 West 120th Street, 

Q i New York, NY 10027-6601 

O 



u 
< 



julesOastro . Columbia. edu 

S. M. Ransom 

Department of Physics, McGill University, Montreal, QC H3A 2T8, Canada 

ABSTRACT 

> 

A Chandra image of the central X-ray source RX J0007. 0+7303 in the supernova 
remnant CTA 1 reveals a point source, a compact nebula, and a bent jet, all of which 
are characteristic of energetic, rotation-powered pulsars. Using the MDM 2.4m tele- 
scope we obtain upper limits in the optical at the position of the point source, (J2000) 
00 h 07 m 01. s 56, +73°03'08."1, determined to an accuracy of 0."1, of B > 25.4, V > 24.8, 
Qlh| R > 25.1, and I > 23.7; these correspond to an X-ray-to-optical flux ratio > 100. 

Neither a VLA image at 1425 MHz, nor a deep pulsar search at 820 MHz using the 
NRAO Green Bank Telescope, reveal a radio pulsar counterpart to an equivalent lu- 
minosity limit at 1400 MHz of 0.02 mJy kpc 2 , which is equal to the lowest luminosity 
known for any radio pulsar. The Chandra point source accounts for ~ 30% of the 
flux of RX J0007. 0+7303, while its compact nebula plus jet comprise « 70%. The 
X-ray spectrum of the point source is fitted with a power-law plus blackbody model 
with r = 1.6 ± 0.6, kToo = 0.13 ± 0.05 keV, and Roq = 0.37 km, values typical of 
a young pulsar. An upper limit of T c °° < 6.6 x 10 5 K to the effective temperature 
of the entire neutron star surface is derived, which is one of the most constraining 
data points on cooling models. The 0.5-10 keV luminosity of RX J0007. 0+7303 is 
~ 4 x 10 31 (d/1.4 kpc) 2 ergs s _1 , but the larger (~ 18' diameter) synchrotron nebula in 
which it is embedded is two orders of magnitude more luminous. These properties allow 
us to estimate, albeit crudely, that the spin-down luminosity of the underlying pulsar 
is in the range 10 36 — 10 37 ergs s" 1 , and support the identification of the high-energy 
7-ray source 3EG J0010+7309 as a pulsar even though its spin parameters have not yet 
been determined. 



Subject headings: ISM: individual (CTA 1) — stars: individual (RX J0007.0+7303, 
3EG J0010+7309) — stars: neutron — supernova remnants 



- 2 - 



1. Introduction 

The supernova remnant CTA 1 (G119.5+10.2) is a radio shell with a diameter of ~ 1?8 (Seward, 
Schmidt, & Slane 1995) and a center-filled X-ray morphology. It has a kinematic distance of 
d = 1.4 ± 0.3 kpc derived from an associated H I shell (Pineault et al. 1993), and the remnant age 
from a Sedov analysis is t$ ~ 1.3 x 10 4 yr (Slane et al. 2004). A detailed X-ray spectral study of 
the SNR using A SCA data showed that the large-scale extended emission is likely of a synchrotron 
origin; its total luminosity is 5.6 x 10 33 (d/1.4 kpc) 2 ergs s _1 in the 0.5-10 keV band (Slane et al. 
1997). The compact ROSAT PSPC source RX J0007.0+7303 is located in the brightest part of the 
synchrotron nebula (Seward et al. 1995; Slane et al. 1997; Brazier et al. 1998). RX J0007.0+7303 
has the requisite properties for a rotation-powered pulsar. Its X-ray luminosity in the 0.1-2.4 keV 
band is 4.3 x 10 31 (d/lA kpc) 2 ergs s _1 (Slane et al. 1997), typical of young pulsars, but atypically 
small compared to the total luminosity of the large synchrotron nebula in CTA 1. Appealing to 
an empirical relation between X-ray luminosity and spin-down power E, Slane et al. (1997) argued 
that the current spin-down power of the pulsar, required to account for the present nonthermal 
emission from the entire synchrotron nebula, is E = 1.7 x 10 36 ergs s _1 . Combining this with 
the Sedov age predicts a rotation period P < 0.17 s and a dipolar surface magnetic field strength 
B p < 6 x 10 12 G. 

Within the boundary of CTA 1 lies 3EG J0010+7309, one of the brighter unidentified EGRET 
sources. With a 95% confidence localization diameter of only 28' and an intermediate Galactic 
latitude of +10?2, this EGRET/SNR coincidence is one of the most convincing among the unproven 
identifications (Brazier et al. 1998), especially because of the likely pulsar in RX J0007. 0+7303. 
The 7-ray source itself has all of the characteristics of the known pulsars. It is not variable; other 
7-ray pulsars show little if any change in flux while most blazars, the other major class of EGRET 
source, are often dramatically variable. Its photon spectral index of 1.85 ± 0.10 (Hartman et al. 
1999) (1.58 ± 0.18 between 70 MeV and 2 GeV, Brazier et al. 1998), is similar to other EGRET 
pulsars, and flatter than that of most blazars. The estimated values of d, ts, E, and P are all typical 
of what one would expect for a young 7-ray pulsar. In particular, the distance to CTA 1 is typical 
for the predicted EGRET pulsar population (e.g., Halpern & Ruderman 1993; Kaaret & Cottam 
1996; Romani & Yadigaroglu 1995; Yadigaroglu & Romani 1995, 1997), and the inferred energetics 
are comparable to those of the Vela pulsar and other sources with similar values of E. 

Recent observations with XMM-Newton resolved RX J0007.0+7303 into a point source and a 
diffuse nebula, and found a two-component (blackbody plus power-law) spectral fit that is consistent 
with a young pulsar (Slane et al. 2004). Here, we present a higher resolution Chandra observation 
of these features that further resolves out a prominent jet-like structure and possible evidence of a 
torus, features typical of pulsar wind nebulae (PWNe) found around the most energetic of pulsars. 
In §2 we present our Chandra imaging and spectroscopic results, in §3 the optical observations, 
and in §4 new radio imaging and deep radio search for pulsations from RX J0007.0+7303. In §5 
we discuss the constraints on E for the putative pulsar in CTA 1 based on its X-ray and (likely) 
7-ray luminosities. Throughout this paper, for clarity, we refer to the small-scale nebula resolved 
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by Chandra within RX J0007.0+7303 as the PWN, as distinct from the large synchrotron nebula, 
which is also likely wind fed. 

2. Chandra Observation of RX J0007.0+7303 

The central region of CTA 1 was observered on 2003 April 13 with the Advanced CCD Imag- 
ing Spectrometer (ACIS; Burke et al. 1997) onboard the Chandra X-ray Observatory (Weisskopf, 
O'Dell, & van Speybroeck 1996). The source RX J0007.0+7303 was positioned at the default lo- 
cation on the back-illuminated S3 chip of the ACIS-S array. The standard TIMED readout with a 
frame time of 3.2 s was used, and the data were collected in VFAINT mode. A total of 50139 s of 
on-time was accumulated, while the effective exposure live-time was 49484 s. All data reduction 
and analysis were performed using the CIAO (V3.0.1), FTOOLS (V5.2), and XSPEC (V11.2) X-ray 
analysis software packages. No time filtering was necessary as the background rate was stable over 
the course of the observation. Photon pile-up was not a consideration as the total count rate in 
the point source of interest was < 0.01 s _1 . 

We used the standard processed and filtered event data with the latest aspect alignments, 
with the exception that the 0.5 pixel (0."25) randomization that is ordinarily applied to the photon 
positions was reversed, restoring slightly sharper images. The 0."5 ACIS pixels slightly undersample 
the on-axis point spread function of the Chandra mirrors in the restored images, as the radius 
encircling 50% of the energy is ~ 0."5. However, the spacecraft dithering and time-dependent 
aspect solution preclude the need for additional randomization, which degrades a study of faint 
point sources embedded in diffuse emission. 

2.1. ACIS Image 

In Figure 1 we compare the Chandra ACIS image of RX J0007.0+7303 from 2003 April 13 
with the combined MOS1 and MOS2 image from the XMM-Newton observation of 2002 February 
21, originally reported by Slane et al. (2004). In the Chandra image, 187 photons (after background 
subtraction) are found within 1" of the source maximum; most are attributable to a point source. 
Figure 2 shows the radial profile of the Chandra source, with a simulated point-spread-function 
(PSF) scaled to the detected counts in the central pixel. There is also a compact surrounding 
nebulosity with 136 photons within a radius of 3" (the PWN), as well as a jet that extends 16" 
to the south, with a « 50° bend to the southwest at 12" from the point source. There is also a 
possible faint extension of the jet for an additional 8" to the southwest, although the latter may be 
part of a general expansion of the jet into a large, low-surface brightness nebulosity extending to 
the west and northwest. The total number of photons in the jet is 45 after background subtraction. 
From the Vela pulsar, similar diffuse emission apparently supplied by its jet has been observed 
(Kargaltsev et al. 2003; Pavlov et al. 2003). 
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In view of Chandra 's unique ability to reveal graphic evidence of a rotation-powered neutron 
star, we refer in this paper to "the pulsar" in RX J0007. 0+7303 rather than the customary cir- 
cumlocution "candidate pulsar." The ACIS frame time of 3.2 s does not permit a search for X-ray 
pulsations in this observation. Prom the XMM-Newton EPIC pn data taken with 6 ms time reso- 
lution, Slane et al. (2004) derive an upper limit of 61% sinusoidal pulsed fraction for frequencies up 
to 83 Hz. Nevertheless, we regard it as appropriate in cases such as this one to classify an X-ray 
source as a pulsar even though its spin parameters are not yet known. The position of the pulsar 
is (J2000) 00 h 07 m 0L s 56, +73 o 03'08."l (see §3 for a description of the astrometry). 

While the similarity between the jets of RX J0007.0+7303 and Vela is striking, it is curious 
that the prominent torus structure of Vela is so much weaker in RX J0007.0+7303, if it is present at 
all. The compact PWN surrounding the pulsar appears to be elongated at an angle perpendicular 
to the inner jet, which by analogy with Vela may plausibly be interpreted as an equatorial torus 
with the jet emerging along the rotation axis. Figure 3 shows the Chandra image in two different 
energy bands, soft (0.2-2.0 keV) and hard (2-8 keV). The PWN is more prominent in the hard 
band, while the pulsar itself is brighter in the soft band, presumably due to the contribution of 
surface thermal emission from the neutron star (see §2.2.2), which has a softer spectrum than the 
synchrotron nebulosity. Apparent spectral differences exist along the length of the jet. The hardest 
regions of the jet are closest to the pulsar, while the softest part of the jet is farther away, below 
the bend. 

The features of the lower-resolution XMM-Newton image are roughly consistent with what one 
would expect from the Chandra image, although the part of the jet past the bend is not obvious 
with XMM-Newton. It is possible that the structure of the jet has changed markedly in the 14 
months between the two observations. Such behavior would be consistent with the rapid changes 
seen in the jet of the Vela pulsar (Kargaltsev et al. 2003; Pavlov et al. 2003), which moves with 
speeds in the range 0.3 — 0.7 c, even in directions transverse to its length. Velocities comparable 
to c were also seen in Chandra observations of Gil. 2-0. 3 (Roberts et al. 2003). At the 1.4 kpc 
distance of CTA 1, the 16" length of the jet corresponds to 3.4 x 10 17 cm, or 130 light days. The 
length of the Vela pulsar's jet is almost the same, « 4.3 x 10 17 (ti/300 pc) cm, which may suggest 
that the two sources have similar spin-down luminosities. The bending and termination of the jet 
can plausibly be explained by its impact against the larger synchrotron nebula. In addition, if the 
pulsar is traveling toward the southeast, the low surface brightness diffuse emission to the west and 
northwest of the pulsar may be due to jet material dumping its energy and dispersing it into the 
larger nebula, as suggested by Pavlov et al. (2003) in the case of Vela. 



2.2. X-ray Spectral Fitting 

Source and background spectra were extracted using the CIAO script psextract to generate 
the photon energy histograms and appropriate response matrices for analysis of three spatial regions: 
the pulsar (r = 1"), the PWN (1" < r < 3"), and the jet (irregular polygon). The PWN region 
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served as background for the pulsar spectrum, while a larger annular background region (11" < 
r < 15") was extracted for use with the PWN and the jet spectra. All spectra were grouped with 
a minimum of 20 counts per bin. 

2.2.1. The PWN and its Intervening Column Density 

We first fitted the 136 background-subtracted counts obtained from the PWN region to an 
absorbed power-law model. This produced a rather unconstrained fit with Nn = 3.5(0.0 — 8.9) x 
10 21 cm" 2 and T = 1.0(0.3 — 1.4), where the la range is given in parentheses. The unabsorbed 
flux in the 0.5 — 10 keV band is 6.6 x 10" 14 ergs cm" 2 s" 1 after adding a presumed contribution, 
scaled by area, from within the pulsar region (r = l"). We note that the best fitted value of 
Nn = 3.5 x 10 21 cm -2 is somewhat larger than expected based on the total Galactic optical 
extinction in the direction of CTA 1. According to the 100 /im IRAS/COBE maps (Schlegel, 
Finkbeiner, & Davis 1998), E(B — V) = 0.414 mag, corresponding to Ay = 1.37, which can be 
considered an upper limit on the extinction to CTA 1. The Predehl & Schmitt (1995) relation 
N H /A V = 1.8 x 10 21 cm" 2 mag" 1 then implies N H = 2.5 x 10 21 cm -2 . The latter value is close 
to that derived by Slane et al. (1997) from joint fits to ROSAT PSPC and ASCA spectra, which 
gave iVn = 2.8(2.3 — 3.4) x 10 21 cm" 2 . Given the distance and intermediate Galactic latitude 
of CTA 1 (+10? 2), it is plausible that the source lies above the bulk of the dust layer, since 
1400sin(10?2) s=s 250 pc. Therefore, we adopt a fixed value of Nu = 2.8 x 10 21 cm" 2 in all further 
spectral fits, the results of which are listed in Table 1. 

The 0.5-10 keV luminosity of the compact PWN region in RX J0007.0+7303 is 1.8x 10 31 ergs s" 1 . 
In comparison, the 0.1-10 keV luminosity of the Vela nebula within a radius of 53" (2.4 x 10 17 cm) 
is (5 - 6) x 10 32 ergs s" 1 (Helfand, Gotthelf, & Halpern 2001; Pavlov et al. 2003), a difference that 
is manifest in the relative prominence of the toroidal features of the Vela nebula. 

2.2.2. The Pulsar 

We next fitted the pulsar spectrum, consisting of 187 background-subtracted counts, to a two- 
component absorbed power-law plus blackbody model, to allow for non-thermal magnetospheric 
emission and thermal emission from the neutron star surface, respectively. The column density was 
held fixed at N n = 2.8 x 10 21 cm" 2 . The best fitted parameters are T = 1.6(1.1 - 2.2) and kT^ = 
0.13(0.08 - 0.18) keV, with unabsorbed fluxes in the 0.5 - 10 keV band of 4.0 x 10" 14 ergs cm" 2 s" 1 
and 8.2 x 10 -15 ergs cm" 2 s" 1 , for the power-law and blackbody components, respectively. The small 
effective radius of the blackbody component, Rqq = 0.37 km, is similar to the value R^ = 0.63 km 
from the VMM-Newton analysis (Slane et al. 2004), and might be evidence of a hot polar cap 
heated by backflow of energetic particles. On the other hand, Slane et al. (2004) also fitted a 
light-element atmosphere model with a larger fixed radius of 10 km, enabling a lower temperature 
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to be derived. Such systematic ambiguities are always present in thermal modeling of neutron star 
spectra. 

Independent of this spectral fit to the pulsar, we also derive an upper limit to the effective 
blackbody temperature of the entire neutron star surface. Assuming d = 1.4 kpc, iVn = 2.8 x 
10 21 cm~ 2 , and a radius at infinity of 12 km, we compared simulated blackbody spectra of increasing 
temperature with the data until the predicted spectrum exceeded the observered counts in the 
lowest energy bin by 3a, or nearly a factor of 2. The resulting upper limit is T c °° < 6.6 x 10 5 K. 
In Figure 4, we compare this value with a range of cooling models that include both standard and 
"exotic" processes (see also Page & Applegate 1992; Tsuruta et al. 2002, and the reviews of neutron 
star cooling by Tsuruta 1998 and Yakovlev & Pethick 2004). Recent calculations show that the 
temperatures or upper limits on all cooling neutron stars, such as Vela, Geminga, 3C 58 (PSR 
J0205+6449), and RX J0007.0+7303, can be fitted with models that include proton and neutron 
superfluidity, and masses < 1.5 M Q (Kaminker, Yakovlev, & Gnedin 2002; Yakovlev et al. 2002), 
as shown in Figure 4. Here, the direct Urea process is operating in the more massive neutron stars. 
Cooling can also be enhanced in a superfluid neutron star by Cooper-pairing neutrino emission, 
which was modeled most recently by Gusakov et al. (2004) and Page et al. (2004). 

2.2.3. The Jet 

For the jet we created a polygonal region which encompassed and followed the contour of 
the jet from 2'.'5 to 22" from the pulsar. A fit to the 45 background-subtracted photons with an 
absorbed power-law model, with Nu again fixed to the ASCA value, yielded V = 1.26(0.86 — 1.69) 
and an unabsorbed flux of 3.1 x 10 -14 ergs cm~ 2 s _1 . While there are probably changes in spectral 
index along the jet that are related to the apparent differences in the soft and hard images (Fig. 
3), there are too few photons to fit spectra to multiple regions. The 0.5-10 keV luminosity of the 
jet, 7.3 x 10 30 ergs s _1 , is almost identical to that of the main (inner plus outer) jet of the Vela 
pulsar, which has Lv(l — 8 keV) ~ 6 x 10 30 ergs s _1 (Pavlov et al. 2003). Within the range of 
uncertainty and variability, these two sources have jets of equal length and luminosity, which is in 
curious contrast to the relative dominance of the toroids in Vela. 



3. Optical Observations of RX J0007.0+7303 

We obtained multiple CCD images of the field of RX J0007.0+7303 using the 2.4m Hiltner 
telescope of the MDM Observatory on 2002 August 26-30. A thinned, back-illuminated SITc 
2048 x 2048 CCD with a spatial scale of 0."275 per 24 ^m pixel was used to cover a 9/4 x 9.'4 
area. The sky conditions were usually photometric enough to calibrate the BVRI images with 
Landolt (1992) standard stars, although moonlight prevented us from reaching the deepest possible 
limiting magnitudes. Total integration times were 100, 72, 120, and 75 minutes in B,V,R, and /, 



- 7- 



respectively. The seeing in the final, summed images, shown in Figure 5, ranged from a best of 0."8 
in / to a worst of l."2 in V. 

We were able to confirm and refine the astrometric accuracy of the Chandra X-ray image 
relative to the optical reference frame of the USNO-A2.0 catalog (Monet et al. 1998) by identifying 
the counterparts of several X-ray sources in the vicinity of RX J0007. 0+7303. X-ray positions and 
counts from the standard data processing, and optical positions and magnitudes from our CCD 
images, are listed in Table 2. The X-ray centroids and source counts are calculated within the 
80% enclosed energy contour of the local point-spread function. All of the optical counterparts 
in Table 2 appear point-like except the first entry, CXOU J000647.2+730109, which is evidently 
a compact galaxy, and perhaps the very faint CXOU J000739. 2+730404, whose structure is not 
clear. The five neighboring sources had mean X-ray-optical offsets of only — 0."ll in right ascension 
and +0."06 in declination, and rms scatter of 0."28 in right ascension and 0."14 in declination. The 
rms scatter is consistent with the statistical accuracies of the individual optical and X-ray source 
positions, while the mean offset is well within the Chandra systematic 90% confidence accuracy of 
0."6. Accordingly, we corrected the Chandra coordinates of the RX J0007.0+7303 point source by 
the small mean offsets, obtaining a final position of (J2000) 00 h 07 m 0L s 56, +73°03'08."1, which we 
expect is accurate to « 0."1 relative to the optical reference frame. 

Limiting (3a) magnitudes at the position of the Chandra point source are B > 25.4, V > 24.8, 
R > 25.1, and / > 23.7. The extinction-corrected limits using the Schlegel et al. (1998) values 
of A B = 1.79, Ay = 1.37, Ar = 1.11, and A/ = 0.80, are B > 23.6, V > 23.4, R > 24.0, and 
/ > 22.9. This corresponds to fx/ fv > 300, or > 83 after correction for extinction. Similar values 
are obtained in the other bands, and are consistent with an isolated neutron star at the distance 
of CTA 1. No extended optical emission is seen associated with the X-ray jet. 

4. Radio Observations of RX J0007.0+7303 
4.1. VLA Image 

The field of RX J0007.0+7303 was imaged at 1425 MHz using the AB configuration of the VLA 
on 2002 May 27, which gave a beam FWHM of 5."9 x 4."4. An exposure time of 1 hour achieved a 
la noise level of 31 ^Jy and a non-detection at the precise Chandra position of the pulsar. This is 
equivalent to a 3a flux density upper limit of 0.09 mJy The relevant part of the image in shown 
in Figure 6. 

4.2. Radio Pulsar Search at GBT 

CTA 1 was searched for pulsations by Nice & Sayer (1997) at radio frequencies of 370 MHz and 
1390 MHz. The sensitivity of those observations, translated to a standard frequency of 1400 MHz 
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using a spectral index a = —1.6 (the median for all pulsars; Lorimer et al. 1995), was ~ 0.2 mJy 
and ~ 0.6 mJy, respectively. Also, Lorimer, Lyne, & Camilo (1998) reached a limit S1400 ~ 0.2 mJy 
in their search at 600 MHz. Our VLA upper limit, S1400 i$ 0.1 mJy, implies a luminosity L1400 = 
Suood 2 < 0.2 mJy kpc 2 . While this is a stringent limit (the least luminous young radio pulsar 
known, in SNR 3C 58, has about twice this luminosity; Camilo et al. 2002c), there are weaker 
pulsars known. We therefore used the extremely sensitive Green Bank Telescope (GBT), where 
CTA 1 is a circumpolar source, to improve on this limit. 

We observed the Chandra position of the pulsar for a total of 19.60 hr on 2003 October 11. We 
used the Berkeley-Caltech Pulsar Machine (BCPM; Backer et al. 1997), an analog/digital filter- 
bank with 96 channels for each of two polarizations, to record a 48 MHz-wide band centered at 
a sky frequency of 820 MHz. After summing in hardware the signals from the two polarizations, 
removing the mean and scaling, the power levels from each frequency channel were recorded with 
4-bit precision to disk every 72 /Lts, for a total of 980 million time samples. [See Camilo et al. (2002c) 
for more details of a similar observation at the GBT using the BCPM.] 

We analyzed the data using standard methods as implemented in the software package PRESTO 
(Ransom 2001). This involved identification and removal of the worst radio- frequency interference 
(much of it generated at the observatory), time-shifting the 96 frequency channels to compensate 
for dispersive interstellar propagation prior to summing the respective samples (doing this once 
for each assumed value of dispersion measure, the integrated column density of free electrons), 
performing a Fast-Fourier Transform of the resulting one-dimensional time series for each trial DM 
after transformation to the solar-system barycenter, and sifting the resulting spectrum for statisti- 
cally significant peaks. The last stage included searching for narrow pulse shapes by summing up 
to eight harmonics of the spectrum; and also searching for moderately accelerated signals, as might 
be expected from a pulsar in a binary system, or a young one with a large period derivative. 

At a distance of 1.4 kpc, both the Taylor & Cordes (1993) and Cordes & Lazio (2002) electron 
density models predict DM ~ 25 cm -3 pc for the putative radio pulsar in CTA 1. However, the 
X-ray-measured Ah ~ 2.8 x 10 21 cm" 2 would predict DM ~ 90 cm -3 pc, assuming an average 
ratio n e /riH ~ 0.1. Also, the maximum Galactic dispersion in the direction of CTA 1 predicted by 
the Cordes & Lazio model is 200 cm" 3 pc, so we searched up to this DM. Although the expected 
period for the pulsar in CTA 1 is P ~ 0.1 s, and almost certainly not less than ~ 20 ms, we chose 
to analyze our data with much higher time resolution in order to retain serendipitous sensitivity to 
any pulsars located within the 15'-diameter GBT beam. In order to simplify the onerous reduction 
task, we first rebinned the data four-fold to a time resolution of 0.288 ms, which we kept for the 
first 240 trial DMs, each separated by 0.5 cm" 3 pc. The effective time resolution at the upper end 
of this range, limited by dispersive smearing within one frequency channel, is about 1 ms. We then 
further halved the time resolution, doubled the DM step, and searched a further 80 time series. 
In the end, we examined « 100 significant signals that appeared at more than one DM, but none 
showed a true interstellar dispersion pattern that would qualify it as a good pulsar candidate. 
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We now estimate the flux density limit of our observation. The ideal rms noise of the observa- 
tion is a = TG^ 1 (nBt)^ 1 / 2 , where T ~ 35 K includes the system temperature on cold sky as well as 
a~7K contribution from Galactic synchrotron emission and ~3K from CTA 1; G = 2 K Jy~ x is 
the telescope gain; n = 2 polarizations; B is the bandwidth; and t is the integration time. With our 
parameters, a « 7 fiJj. Including a ~ 20% factor accounting for losses due to hardware limitations 
and a threshold signal-to-noise ratio of 8, we obtain an estimated sensitivity for a sinusoidal pulse 
shape of ~ 65 ^uJy For a pulse with duty cycle to = 0.1 P, this limit is further improved by a factor 
of ~ 3. Translating this to 1400 MHz with a ~ —1.6, we obtain S1400 < 30 ^Jy for a sinusoidal 
pulse shape or 51400 ^5 10 /"Jy f° r w = 10%. Strictly, these limits apply to relatively long pulse 
periods; for P < 50 ms the limits are somewhat worse, the more so as DM/P increases. 

Taking 10 /iJy as the flux density limit at 1400 MHz for the pulsar in CTA 1 implies L1400 ^ 0.02 
mJy kpc 2 . This is equal to the lowest luminosity known for a radio pulsar, and is a factor of ss 20 
below the least-luminous pulsar known with an age of less than 1 Myr. If the pulsar in CTA 1 
emits radio waves that intersect the Earth, its luminosity is lower than that of any detected radio 
pulsar; perhaps more likely, if it is a radio pulsar at all, its beam does not intersect the Earth. 

5. Discussion: What is the Pulsar El 

The most important unknown property of RX J0007.0+7303 is its spin-down luminosity. Pre- 
dictions of E typically make use of correlations between Lx and E among known pulsars (e.g., 
Seward & Wang 1988; Possenti et al. 2002), but this method is fraught with uncertainty, as well as 
special difficulties in the case of RX J0007.0+7303. The 0.5-10 keV luminosity of RX J0007.0+7303 
that we measure with Chandra, including the point and diffuse components, is 4.0 x 10 31 ergs s _1 . 
This is slightly smaller than the value of 5.4 x 10 31 ergs s _1 measured by Slane et al. (2004) using 
XMM-Newton, the difference probably attributable to an additional contribution of diffuse emis- 
sion contained in the larger XMM-Newton extraction region. The spin-down luminosity predicted 
by Slane et al. (2004) using the Seward & Wang (1988) relation is then E ~ 5 x 10 34 ergs s -1 , com- 
parable to the values found for older 7-ray pulsars such as Geminga, which has an E = 3.3 x 10 34 
ergs s -1 . This interpretation is troubling in its neglect of the large synchrotron nebula, for which 
there is no proposed power source other than the pulsar. It is impossible to measure the entire syn- 
chrotron nebula of CTA 1 with Chandra because of its large size, low surface brightness, and lack 
of a clear boundary. But since the total synchrotron emission seen by ASCA exceeds the compact 
source luminosity by more than two orders of magnitude, consideration of the total energy budget 
would favor a minimum value of E ~ 1.7 x 10 36 ergs s _1 , as originally proposed by Slane et al. 
(1997). 

We may also look to more detailed X-ray studies of pulsars with PWNe to estimate the spin- 
down power of CTA 1. From an empirical perspective, bright PWNe are found only around pulsars 
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with E > 3 x 10 36 ergs s _1 (Gotthelf 2004). 1 Above this threshold, PWN emission typically 
exceeds the pulsar point-source luminosity by an order of magnitude. Furthermore, the spectral 
index of the pulsar in the 2-10 keV band correlates with E according to the Gotthelf (2003) relation 
r PSR = 2.08- 0.29 E~z /2 , where E 38 is the spin-down power in units of 10 38 ergs s 1 Applying this 
to CTA 1 for which Tpsr = (1.1 — 2.2) suggests that E > 6 x 10 36 ergs s" 1 , allowing for scatter in 
the relation. 

Further indication that E is in this higher range is suggested by the 7-ray luminosity of 
3EG J0010+7309, which provides a lower-limit on the spin-down power. At d = 1.4 kpc, the 7-ray 
flux integrated above 100 MeV corresponds to L 7 = 6 x 10 34 ergs s _1 if isotropic, or 5 x 10 33 
ergs s _1 if beamed into 1 steradian as is often (but arbitrarily) assumed. Thus, a prediction of 
E ~ 5 x 10 34 ergs s _1 is incompatible with isotropic 7-ray emission, and requires a substantial 
degree of beaming. On the other hand, if E ~ 6 x 10 36 ergs s _1 , then the ratio L^/E is only 0.01 in 
the isotropic case, which is similar to Vela, whose E = 6.9 x 10 36 ergs s _1 . Note that this argument 
is independent of whether the X-ray spectrum of RX J0007. 0+7303 can be extrapolated over four 
orders of magnitude in energy to join the 7-ray spectrum. Such a connection was suggested by 
Slane et al. (2004), but it is impossible to verify, as well as unnecessary for associating the sources. 
The 0.5 — 10 keV power law of the Geminga pulsar clearly does not connect with its 7-ray spectrum 
(Jackson et al. 2002); rather, they are separate spectral components arising from different emission 
mechanisms. 



6. Conclusions 

The resemblance of the X-ray morphology of RX J0007.0+7303 to the PWN and jets seen 
in Chandra images of the Crab (Weisskopf et al. 2000), Vela (Helfand et al. 2001; Pavlov et al. 
2001, 2003), PSR B1509-58 (Gaensler et al. 2002), and other newly discovered young pulsars 
(Halpern et al. 2001; Lu et al. 2002; Camilo et al. 2002a; Hessels et al. 2004), enables us to point 
confidently to "the pulsar" in the SNR CTA 1 even though its spin parameters have not yet been 
determined. Similar conclusions were drawn from Chandra observations of supernova remnants 
GO.9+0.1 (Gaensler, Pivovaroff, & Garmire 2001), IC 443 (Olbert et al. 2001), 3C 396 (Olbert et 
al. 2003), and another unidentified EGRET source, GeV J1809-2327 (Braje et al. 2002). Using 
clues gathered from its X-ray luminosity, morphology, spectrum, and likely identification with 
3EG J0010+7309, we predict that the spin-down power of the pulsar in RX J0007.0+7303 is in the 
range 10 36 — 10 37 ergs s _1 . 

The upper limit of T c °° < 6.6 x 10 5 K that we derive on the effective temperature of the full 
neutron star surface in CTA 1 is more constraining of cooling models than nearly all other cases. 



1 Extended emission recently found around the nearby (d ~ 160 pc) Geminga pulsar (Caraveo et al. 2003) has Lx 
of only ~ 6.5 x 10 28 ergs s _1 , a mere 2 x 10~ 6 of its E or 2% of its Lx- Such weak emission would not have been 
detected at the distance of CTA 1. 
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As described by Yakovlev et al. (2002), X-ray observations do not yet require exotic phases of 
matter such as pion and kaon condensates, or free quarks. However, it is necessary to allow a 
range of neutron star masses if the measured temperatures and ages are to be accommodated by 
a single equation of state and theory of superfluid properties (Gusakov et al. 2004). In this sense, 
the neutron star in CTA 1 is cooler for its age than most others, and similar to Vela and Geminga 
in requiring a mass > 1.42 M . If so, it may become an issue to understand why single neutron 
stars often have inferred masses greater than the precise values that are measured in binary radio 
pulsars, which cluster tightly around 1.35 M© (Thorsett k Chakrabarty 1999; Stairs et al. 2002; 
Lyne et al. 2004; Weisberg k Taylor 2003). Discussion of that potential problem is beyond the 
scope of this paper. 

Energetic pulsars are being discovered with extremely small radio luminosities; whether this is 
intrinsic or due to beaming remains unclear. Recent examples include PSR J2229+6114 coincident 
with 3EG J2227+6122 (Halpern et al. 2001), PSR J1930+1852 in G54.1+0.3 (Lu et al. 2002; 
Camilo et al. 2002a), PSR J1124-5916 in G292.0+1.8 (Hughes et al. 2001; Camilo et al. 2002b), 
and PSR J0205+6449 in 3C 58 (Murray et al. 2002; Camilo et al. 2002c). All of these pulsars have 
E = (1 — 3) x 10 37 ergs s _1 and luminosities at 1400 MHz of ~ 1 mJy kpc 2 . Such radio luminosities 
are less than those of ~ 95% of all pulsars, and less than every other pulsar younger than 2 x 10 yr 
(Camilo et al. 2002b). In comparison, the luminosity upper limit for RX J0007. 0+7303 is only 
k, 0.02 mJy kpc 2 , which makes it fainter in the radio than all pulsars except the nearby Geminga 
and perhaps the putative pulsar powering 3EG J1835+5918 (Mirabal et al. 2000; Mirabal k Halpern 
2001; Halpern et al. 2002). If due to unfavorable orientation with respect to the radio beam, the 
likely existence of a radio-quiet 7-ray pulsar in CTA 1 tends to support the outer-gap model, in 
which the 7-rays arise from regions far out in the magnetosphere and are emitted into a wide fan 
beam, illuminating a much larger fraction of the sky than does the narrow radio beam. 

We thank Don Backer for making the BCPM available to the GBT user community. This is 
work was supported by SAO grant GO3-4064X. 
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Fig. 1 — XMM-Newton and Chandra images of the central source RX J0007.0+7303 in CTA 1 in 
the 0.5 — 8 keV band, at two different display scales. North is up and east is to the left. Top panels 
are 256" x 256". Bottom panels are 64" x 64". Left: Combined XMM-Newton MOS1 and MOS2 
images, smoothed with a Gaussian of a = 1". Right: Chandra ACIS-S3 image, smoothed with a 
Gaussian of a = 0."5. The position of the point source is (J2000) 00 h 07 m 0L s 56, +73°03'08"1. [A 
color version of this figure is available from the author.] 




Fig. 2. — Chandra radial profile of RX J0007. 0+7303 (histogram) compared to the PSF (dashed 
line), which is calculated for the observed spectral distribution. The PSF is scaled to match the 
detected counts in the central pixel and the background measured in the annulus 25" < r < 30". 
A clear excess is seen corresponding to the faint nebulosity at r < 3", and additional enhancement 
at 3" < r < 20" is likely due to the jet and fainter diffuse emission. 



Fig. 3. — Chandra images of the central source RX J0007.0+7303 in CTA 1 in different energy 
bands and smoothed with a Gaussian of a = 0."5. North is up and east is to the left. Left: Soft 
(0.2-2.0 keV). Right Hard (2-8 keV). Top panels are 64" x 64". Bottom panels are 32" x 32". [A 
color version of this figure is available from the author.] 
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Fig. 4. — Left: Upper limit on the surface temperature of the neutron star in CTA1, in comparison 
with cooling models of Page (1998). The horizontal line denotes a possible age range of (1 — 3) x 
10 4 yr. This figure is reproduced from Slane, Helfand, & Murray (2002), and includes their upper 
limit on 3C 58. Right: This figure is taken from Yakovlev et al. (2002), and includes cooling models 
calculated assuming strong proton superfluidity and weak neutron superfluidity. Here, the direct 
Urea process is not completely suppressed, and operates in the more massive neutron stars. 
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Fig. 5.— CCD images at the location of the central X-ray source RX J0007.0+7303 in CTA 1 
obtained with the 2.4m Hiltner telescope. North is up, and east is to the left. Each panel is 70" 
on a side. Seeing ranges from 0."8 in / to l."2 in V. Interference fringes in the /-band image 
are artifacts. Limiting (3<r) magnitudes at the position of the Chandra point source, (J2000) 
00 h 07 m 01. s 56, +73°03'08."l as indicated by the crosses, are B > 25.4, V > 24.8, R > 25.1, and 
/ > 23.7. [A higher resolution version of this figure is available from the author.] 
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QQ^ 07 m 30 s 15 s Q 7 m 0Q s 45 s 3Q s 15 s 

R. A. (J2000) 

Fig. 6. — VLA image of the central region of CTA 1 at 1425 MHz, taken in AB configuration 
with a beam FWHM of 5."9 x 4."4. The cross marks the Chandra location of the neutron star in 
RX J0007. 0+7303. The la noise level is 31 fiJy per beam at this position. For comparison, the 
point source marked "S" has a flux density of 180 /xJy. The large-scale periodic noise in the image 
is an artifact. 
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Table 1. Spectral Fits to Chandra Observation of RX J0007.0+7303 



Parameter 51 Pulsar PWN Jet 

Net counts 187 136 45 

N K (cm~ 2 ) (fixed) 2.8 x 10 21 2.8 x 10 21 2.8 x 10 21 

T 1.62(1.10 - 2.20) 0.97(0.61 - 1.35) 1.26(0.86 - 1.69) 

C p i (photons cm" 2 s" 1 keV" 1 ) 13 5.8 x 10~ 6 4.3 x 10~ 6 2.9 x lO" 6 

Fpi (ergs cm" 2 s~ 1 ) c 4.0 x 10" 14 7.5 x 10~ 14 3.1 x 10" 14 

L p i (ergs cm" 2 s" 1 ) c 9.5 x 10 30 1.8 x 10 31 7.3 x 10 30 

/cToo (keV) d 0.13(0.08 - 0.18) 

iioo (km) d 0.37 

F bb (ergs cm" 2 s^ 1 ) c 8.2 x 10" 15 

Loo(bol) (ergs cm" 2 s- 4 ) d 5.0 x 10 30 

X 2 [dof] 0.58[6] 0.68[4] 0.75 [5] 

a Uncertainties are 68% confidence for one interesting parameter. 
b Power-law normalization at 1 keV. 

c Unabsorbed fluxes and luminosities in the 0.5-10 keV band, assuming d = 1.4 kpc. 
d Spherical blackbody parameters for d = 1.4 kpc. 



Table 2. Chandra Sources and Optical Counterparts 



Source 


X-ray Position 


Photons 


Optical Position 


A R.A. 


A Dccl. 


B 


R 




R.A. Decl. 


(0.1-10 keV) 


R.A. Dccl. 


(") 


(") 


(mag) 


(mag) 


CXOU J000647.2+730109 


00 06 47.244 +73 01 09.12 


15 


00 06 47.369 +73 01 09.33 


-0.55 


-0.21 


24.8 


22.7 


CXOU J000659. 4+730219 


00 06 59.414 +73 02 19.16 


23 


00 06 59.393 +73 02 19.00 


+0.09 


+0.16 


20.3 


18.4 


RX J0007.0+7303 


00 07 01.537 +73 03 08.17 


274 








> 25.4 


> 25.1 


CXOU J000739.2+730404 


00 07 39.276 +73 04 04.74 


36 


00 07 39.289 +73 04 04.58 


-0.06 


+0.16 


25.2 


23.6 


CXOU J000740.3+730649 


00 07 40.343 +73 06 49.43 


176 


00 07 40.408 +73 06 49.31 


-0.28 


+0.12 


23.3 


21.6 


CXOU J000741. 9+730227 


00 07 41.936 +73 02 27.17 


158 


00 07 41.877 +73 02 27.11 


+0.26 


+0.06 


23.6 


22.1 



Note. - 



Units of right ascension are hours, minutes, and seconds, and units of declination are degrees, arcminutcs, and arcseconds. 



